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Abstract 

Using 34.7 pb^^ of data collected with the LHCb detector, the inclusive production 
of the X(3872) meson in pp collisions at ^/s = 7 TeV is observed for the hrst time. 
Candidates are selected in the X(3872) — )• J / ipir'^ it~ decay mode, and used to 
measure 

a{pp X(3872) + anything) S(X(3872) ^ J/^vr+vr") = 

5.4 ± 1.3 (stat) ± 0.8 (syst) nb , 

where a{pp — )■ X(3872) + anything) is the inclusive production cross-section of 
X(3872) mesons with rapidity in the range 2.5 — 4.5 and transverse momentum in 
the range 5 — 20GeV/c. In addition the masses of both the X(3872) and 'ip{2S) 
mesons, reconstructed in the J/'07r+7r~ final state, are measured to be 

"ix(3872) = 3871.95 ± 0.48 (stat) ± 0.12 (syst) MeV/c^ and 
m^(25) = 3686.12 ±0.06 (stat) ±0.10 (syst) MeV/c^ 
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1 Introduction 



The X(3872) particle was discovered in 2003 by the Belle collaboration in the — )■ 
X(3872)ir±, X(3872) ^ J/V'TT+Tr" decay chain [I]. Its existence was confirmed by the 
CDF p], D0 mj and BaBar [1] collaborations. The discovery of the X(3872) particle and 
the subsequent observation of several other new states in the mass range 3.9 — 4.7 GeV/ 
have led to a resurgence of interest in exotic meson spectroscopy [5j. 

Several properties of the X(3872) have been determined, in particular its mass [HIITIIH] 
and the dipion mass spectrum in the decay X(3872) — )■ J/ip7i~^ii~ [9l [7j, but its quantum 
numbers, which have been constrained to be either J^'-^ = 2 ^ or l"*""*" [10], are still 
not established. Despite a large experimental effort, the nature of this new state is still 
uncertain and several models have been proposed to describe it. The X(3872) could 
be a conventional charmonium state, with one candidate being the 7]c2{^D) meson [5]. 
However, the mass of this state is predicted to be far below the observed X(3872) mass. 
Given the proximity of the X(3872) mass to the D*^D^ threshold, another possibility is 
that the X(3872) is a loosely bound D*^D^ 'molecule', i.e. a ((mc)(cm)) system For 
this interpretation to be valid the mass of the X(3872) should be less than the sum of D*^ 
and masses. A further, more exotic, possibility is that the X(3872) is a tetraquark 
state [Hj. 

Measurements of X(3872) production at hadron colliders, where most of the produc- 
tion is prompt rather than from 6-hadron decays, may shed light on the nature of this 
particle. In particular, it has been discussed whether or not the possible molecular nature 
of the X(3872) is compatible with the production rate observed at the Tevatron [121 US]- 
Predictions for X(3872) production at the LHC have also been published [13j. 

This paper reports an observation of X(3872) production in pp collisions at ^/s = 
7 TeV using an integrated luminosity of 34.7 pb~^ collected by the LHCb experiment. The 
X(3872) —7- J/tpTT^TT' selection is optimized on the similar but more abundant ip{2S) — t- 
J/ip'7i~^7i~ decay. The observed X(3872) signal is used to measure both the X(3872) 
mass and the production rate from all sources including 6-hadron decays, i.e. the absolute 
inclusive X(3872) production cross-section in the detector acceptance multiplied by the 
X(3872) — )■ J/ilJ7r~^n~ branching fraction. 

2 The LHCb spectrometer and data sample 

The LHCb detector is a forward spectrometer [13] at the Large Hadron Collider (LHC). It 
provides reconstruction of charged particles in the pseudorapidity range 2 < i] < 5. The 
detector elements are placed along the LHC beam line starting with the vertex detector 
(VELO), a silicon strip device that surrounds the proton-proton interaction region. It 
is used to reconstruct both the interaction vertices and the decay vertices of long-lived 
hadrons. It also contributes to the measurement of track momenta, along with a large 
area silicon strip detector located upstream of a dipole magnet and a combination of sili- 
con strip detectors and straw drift-tubes placed downstream. The magnet has a bending 
power of about 4 Tm. The combined tracking system has a momentum resolution 6p/p 
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that varies from 0.4% at 5 GeV/c to 0.6% at 100 GeV/c. Two ring imaging Cherenkov 
(RICH) detectors are used to identify charged hadrons. The detector is completed by elec- 
tromagnetic calorimeters for photon and electron identification, a hadron calorimeter, and 
a muon system consisting of alternating layers of iron and multi-wire proportional cham- 
bers. The trigger consists of a hardware stage, based on information from the calorimeter 
and muon systems, followed by a software stage which applies a full event reconstruction. 

The cross-section analysis described in this paper is based on a data sample collected 
in 2010, exclusively using events that passed dedicated J/ip trigger algorithms. These 
algorithms selected a pair of oppositely-charged muon candidates, where either one of the 
muons had a transverse momentum larger than 1.8 GeV/c or one of the two muons 
had pt > 0.56 GeV/c and the other > 0.48 GeV/c. The pair of muons was required 
to originate from a common vertex and have an invariant mass in a wide window around 
the J/tp mass. The X(3872) mass measurement also uses events triggered with other 
algorithms, such as single-muon triggers. To avoid domination of the trigger CPU time 
by a few events with high occupancy, a set of cuts was applied on the hit multiplicity 
of each sub-detector used by the pattern recognition algorithms. These cuts reject high- 
multiplicity events with a large number of pp interactions. 

The accuracy of the X(3872) mass measurement relies on the calibration of the track- 
ing system [15] . The spatial alignment of the tracking detectors, as well as the calibration 
of the momentum scale, are based on the J/^/^ — )■ fi~^fi~ mass peak. This was carried 
out in seven time periods corresponding to known changes in the detector running con- 
ditions. The procedure takes into account the effects of QED radiative corrections which 
are important in this decay. 

The analysis uses fully simulated samples based on the Pythia 6.4 generator [16] 
configured with the parameters detailed in Ref. ^7\. The EvtGen [18], Photos [19] 
and Geant4 [20] packages are used to describe the decays of unstable particles, model 
QED radiative corrections and simulate interactions in the detector, respectively. The 
X(3872) — i- J/ipTT^TT^ Monte Carlo events are generated assuming that the p resonance 
dominates the dipion mass spectrum, as established by the CDF [9J and Belle [7J data. 

3 Event selection 

To isolate the X(3872) signal, tight cuts are needed to reduce combinatorial background 
where a correctly reconstructed J/ip meson is combined with a random tt+tt" pair from 
the primary pp interaction. The cuts are defined using reconstructed ip{2S) — ?■ J/ip-K^-K' 
decays, as well as 'same-sign pion' candidates satisfying the same criteria as used for the 
X(3872) and ip{2S) selection but where the two pions have the same electric charge. The 
KuUback-Leibler (KL) distance [21] is used to suppress duplicated particles created by 
the reconstruction: if two particles have a symmetrized KL divergence less than 5000, 
only that with the higher track fit quality is considered. 

J/ip ^ yU+yU^ candidates are formed from pairs of oppositely-charged particles identi- 
fied as muons, originating from a common vertex with a per degree of freedom (x^/ndf) 



2 



smaller than 20, and with an invariant mass in the range 3.04 — 3.14 GeV/c^. The two 
muons are each required to have a momentum above 10 GeV/c and a transverse momen- 
tum above 1 GeV/c. To reduce background from the decay in flight of pions and kaons, 
each muon candidate is required to have a track fit x^/ndf less than 4. Finally J/ip 
candidates are required to have a transverse momentum larger than 3.5 GeV/c. 

Pairs of oppositely-charged pions are combined with J/ip candidates to build ip{2S) 
and X(3872) candidates. To reduce the combinatorial background, each pion candidate 
is required to have a transverse momentum above 0.5 GeV/c and a track fit x^/ndf less 
than 4. In addition, kaons are removed using the RICH information by requiring the 
likelihood for the kaon hypothesis to be smaller than that for the pion hypothesis. A 
vertex fit is performed p2| that constrains the four daughter particles to originate from a 
common point and the mass of the muon pair to the nominal J / if) mass |23] ■ This fit both 
improves the mass resolution and reduces the sensitivity of the result to the momentum 
scale calibration. To further reduce the combinatorial background the x^/ndf of this fit 
is required to be less than 5. Finally, the requirement Q < 300 MeV/c^ is applied where 
Q = M^^TJT, — M^^ — Mtt^, and M^^^^, M^^ and M^^^ are the reconstructed masses before 
any mass constraint; this requirement removes 35% of the background whilst retaining 
97% of the X(3872) signal. 

Figure [l] shows the J/?/'7r+7r~ mass distribution for the selected candidates, with clear 
signals for both the ip{2S) and the X(3872) mesons, as well as the J/ipTc^iT^ mass distri- 
bution of the same-sign pion candidates. 



4 Mass measurements 

The masses of the ip{2S) and X(3872) mesons are determined from an extended un- 
binned maximum likelihood fit of the reconstructed J/^/'tt'^tt" mass in the interval 
3.60 < Mj/^^TT < 3.95 GeV/c^. The ip{2S) and X(3872) signals are each described with a 
non-relativistic Breit-Wigner function convolved with a Gaussian resolution function. The 
intrinsic width of the ^(2^) is fixed to the PDG value, r^(2s) = 0.304 MeV/c^ [23]. The 
Belle collaboration recently reported [7] that the X(3872) width is less than 1.2 MeV/c^ 
at 90% confidence level; we fix the X(3872) width to zero in the nominal fit. The ratio of 
the mass resolutions for the X(3872) and the ip{2S) is fixed to the value estimated from 
the simulation, o"x(3872)/'^^{25) ~ ^■<ii- 

Studies using the same-sign pion candidates show that the background shape can 
be described by the functional form /(M) oc (M — mthY° exp{—ciM — C2M^), where 
mth = mj/^ + 2m,r = 3376.05 MeV/c^ [23J is the mass threshold and cq, ci and C2 are 
shape parameters. To improve the stability of the fit, the parameter C2 is fixed to the 
value obtained from the same-sign pion sample. 

In total, the fit has eight free parameters: three yields {ip{2S), X(3872) and back- 
ground), two masses {ip{2S) and X(3872)), one resolution parameter, and two background 
shape parameters. The correctness of the fitting procedure has been checked with simpli- 
fied Monte Carlo samples, fully simulated Monte Carlo samples, and samples containing 
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Figure 1: Invariant mass distribution of J/^vr^vr^ (points with statistical error bars) and same- 
sign J /ipTT^TT^ (filled histogram) candidates. The curves are the result of the fit described in 
the text. The inset shows a zoom of the X(3872) region. 



Table 1: Results of the fit to the J /ip-K^Ti invariant mass distribution of Fig. [T 



Fit parameter or derived quantity 


^{2S) 


X(3872) 


Number of signal events 


3998 ± 83 




565 ± 62 


Mass m [ MeV/c^] 


3686.10 ±0.06 


3J 


^71.88±0.48 


Resolution a [ MeV/c^] 


2.54 ±0.06 




3.33 ±0.08 


Signal-to-noise ratio in ±3(T window 


1.5 




0.15 


Number of background events 


73094 d 


z2l 


^2 



a mixture of fully simulated Monte Carlo signal events and same-sign background events 
taken from the data. The fit results are shown in Fig. [T] and Table [TJ The fit does not 
account for QED radiative corrections and hence underestimates the masses. Using a 
simulation based on Photos [IH] the biases on the X(3872) and ip{2S) masses are found 
to be -0.07 ± 0.02 MeV/c^ and -0.02 ± 0.02 MeV/c^, respectively. The fitted mass values 
are corrected for these biases and the uncertainties propagated in the estimate of the 
systematic error. 

Several other sources of systematic effects on the mass measurements are considered. 
For each source, the complete analysis is repeated (including the track fit and the momen- 
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turn scale calibration when needed) under an alternative assumption, and the observed 
change in the central value of the fitted masses relative to the nominal results assigned 
as a systematic uncertainty. The dominant source of uncertainty is due to the calibration 
of the momentum scale. Based on checks performed with reconstructed signals of various 
mesons decaying into two-body final states (such as tt+tt", K^n'^ and fJ,~^n~) a relative 
systematic uncertainty of 0.02% is assigned to the momentum scale which translates 
into a 0.10 (0.08) MeV/c^ uncertainty on the X(3872) (^^(25*)) mass. After the calibration 
procedure with the J/ip ^ Z^"*"/^" decay, a ±0.07% variation of the momentum scale re- 
mains as a function of the particle pseudorapidity r/. To first order this effect averages out 
in the mass determination. The residual impact of this variation is evaluated by parame- 
terizing the momentum scale as function of r] and repeating the analysis. The systematic 
uncertainty associated with the momentum calibration indirectly takes into account any 
effect related to the imperfect alignment of the tracking stations. However, the alignment 
of the VELO may affect the mass measurements through the determination of the hori- 
zontal and vertical slopes of the tracks. This is investigated by changing the track slopes 
by amounts corresponding to the 0.1% relative precision with which the length scale along 
the beam axis is known . Other small uncertainties arise due to the limited knowledge 
of the X(3872) width and the modelling of the resolution. The former is estimated by 
fixing the X(3872) width to 0.7 MeV/c^ instead of zero, as suggested by the likelihood 
published by Belle [7]. The latter is estimated by fixing the ratio (Jx{3872)/<^ip{2S) using 
the covariance estimates returned by the track fit algorithm on signal events in the data 
sample, rather than using the mass resolutions from the simulation. The effect of back- 
ground modelling is estimated by performing the fit on two large samples, one with only 
Monte Carlo signal events, and one containing a mixture of Monte Carlo signal events 
and background candidates obtained by combining a. J/ip candidate and a same-sign pion 
pair from different data events: the difference in the fitted mass values is taken as a 
systematic uncertainty. The amount of material traversed in the tracking system by a 
particle is estimated to be known to a 10% accuracy [25]; the magnitude of the energy 
loss correction in the reconstruction is therefore varied by 10%. The assigned systematic 
uncertainties are summarized in Table [2] and combined in quadrature. 

Systematic checks of the stability of the measured ip{2S) mass are performed, splitting 
the data sample according to different run periods or to the dipole magnet polarity, or 
ignoring the hits from the tracking station before the magnet. In addition, the measure- 
ment is repeated in bins of the p, and Q values of the ip{2S) signal. No evidence for a 
systematic bias is found. 

5 Determination of the production cross-section 

The observed X(3872) signal is used to measure the product of the inclusive production 
cross-section a{pp — )• X(3872) + anything) and the branching fraction i3(X(3872) — )• 
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Table 2: Systematic uncertainties on the ip{2S) and X(3872) mass measurements. 



Category 


Source of uncertainty 


Am [ 
^Pi2S) 


MeV/c2] 
X(3872) 


Mass fitting 


AT i- 1 * Jxl 

iNlatural widtfi 
Radiative tail 
Resolution 
Background model 


0.02 
0.02 


0.01 
0.02 
0.01 
0.02 


Momentum calibration 


Average momentum scale 

f] dependence of momentum scale 


0.08 
0.02 


0.10 
0.03 


Detector description 


Energy loss correction 


0.05 


0.05 


Detector alignment 


Track slopes 


0.01 


0.01 


Total 




0.10 


0.12 



J/ip7i~^Tr ), according to 

a(pp ^ A-(3872) + auythmg) B(A-(3872) ^ J/y,.-.-) = ^^(J/^ )£,... ' 

where N^^^^^^-^ is the efficiency-corrected signal yield, ^ is a correction factor to the 
simulation-derived efficiency that accounts for known differences between data and sim- 
ulation, B{J/ijj — )■ /i"*"//") = (5.93 ± 0.06)% [23j is the J/ijj — )■ branching fraction, 
and £int is the integrated luminosity. 

The absolute luminosity scale was measured at specific periods during the 2010 data 
taking |26] using both Van der Meer scans [27] and a beam-gas imaging method [28] . The 
instantaneous luminosity determination is then based on a continuous recording of the 
multiplicity of tracks in the VELO, which has been normalized to the absolute luminosity 
scale f2B] . The integrated luminosity of the sample used in this analysis is determined to 
be £int = 34.7 ± 1.2 pb~^, with an uncertainty dominated by the knowledge of the beam 
currents. 

Only X(3872) candidates for which the J/ip triggered the event are considered, keeping 
70% of the raw signal yield used for the mass measurement. In addition, the candidates 
are required to lie inside the fiducial region for the measurement, 

2.5 < y < 4.5 and 5 < < 20 GeV/c , (2) 

where y and are the rapidity and transverse momentum of the X(3872). This region 
provides a good balance between a high efficiency (92% of the triggered events) and a low 
systematic uncertainty on the acceptance correction. 

The corrected yield A^x°(3872) ~ 9140 ± 2224 is obtained from a mass fit in the narrow 
region 3820 — 3950 MeV/c^, with a linear background model and the same X(3872) signal 
model as used previously but with the mass and resolution fixed to the central values 
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Table 3: Relative systematic uncertainties on the X(3872) production cross-section measure- 
ment. The total uncertainty is the quadratic sum of the individual contributions. 



OUUiUc Ui UliCci LdliiLj 


Arr/rr \%] 
/-\u j u [ /oj 


^(3872) polarization 


1 


V ( 'i'>i70\ Acrritr mr^AcA 

A [oo (A) decay moflei 


1 n 
i.U 


yvyooizj uecay wiuin 


o.u 


ividoo 1 COVJI U. LlUii 




Background model 


6.4 


Tracking efficiency 


7.4 


Track cut 


2.0 


Vertex cut 


3.0 


Muon trigger efficiency 


2.9 


Hit-multiplicity cuts 


3.0 


Muon identification 


1.1 


Pion identification 


4.9 


Integrated luminosity 


3.5 


J/i/j — )■ fi'^n~ branching fraction 


1.0 


Total 


14.2 



presented in Sect. |4j In this ffi, each candidate is given a weight equal to the reciprocal 
of the total signal efficiency estimated from simulation for the y and of that candi- 
date. A second method based on the sWeight [21] technique was found to give consistent 
results. The average total signal efficiency in the fiducial region of Eq. [2] is estimated to 
be ^x{?,Ki2)l^xjl,-^i2) ~ 4-2%, where A^x(3872) is the observed signal yield obtained from a 
mass ffi without weighting the events. This low value of the efficiency is driven by the 
geometrical acceptance and the requirement on the of the J j meson. 

The quantity ^ of Eq. [l]is the product of three factors. The ffist two, 1.024±0.011 [30] 
and 0.869 ± 0.043, account for differences between the data and simulation for the effi- 
ciency of the muon and pion identifications, respectively. The third factor, 0.92 ± 0.03, 
corresponds to the efficiency of the hit-multiplicity cuts applied in the trigger, which is 
not accounted for in the simulation. It is obtained from a ffi of the distribution of the 
number of hits in the VELO. 

The relative systematic uncertainties assigned to the cross-section measurement are 
listed in Table [3| and quadratically add up to 14.2%. The cross-section measurement is 
performed under the most favoured assumption for the quantum numbers of the X(3872) 
particle, J^^ = l'^^ [31], which is used for the generation of Monte Carlo events. No 
systematic uncertainty is assigned to cover other cases. Besides the uncertainties already 
mentioned on B{J/ip — ^ f^~^f^~)y '^mt and ^, the following sources of systematics on iV'x(3872) 
are considered. The dominant uncertainty is due to differences in the efficiency of track 
reconstruction between the data and simulation. This is estimated to be 7.4% using a 
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data driven tag and probe approach based on J/ip /^^/^^ candidates. An additional 
uncertainty of 0.5% per track is assigned to cover differences in tlie efficiency of the track 
X^/ndf cut between data and simulation. Similarly, a 3% uncertainty is assigned due to 
the effect of the vertex cuts. 

Other important sources of uncertainty are due to the modelling of the signal and 
background mass distributions. Repeating the mass fit with the X(3872) decay width 
fixed to 0.7 MeV/c^ instead of zero results in a 5% change of the signal yield. Similarly, 
the uncertainties due to the X(3872) mass resolution are estimated by repeating the mass 
fit with different fixed mass resolutions: first changing it by the statistical uncertainty 
reported in Table [T| and then changing it by the systematic uncertainty resulting from the 
knowledge of the resolution ratio crx(3872)/o'^{2S)^ as described in Sect. |4j The combined 
effect on the X(3872) signal yield corresponds to a 2.5% systematic uncertainty. 

Using an exponential rather than linear function to describe the background leads to 
a change of 6.4% in signal yield, which is taken as an additional systematic uncertainty. 

The unknown X(3872) polarization affects the total efficiency, mainly through the J/i/j 
reconstruction efficiency. The dipion system is less affected, in particular the efficiency 
is found to be constant as a function of the dipion mass. The simulation efficiency, 
determined assuming no J/ip polarization, is recomputed in two extreme schemes for the 
J/^p polarization (fully transverse and fully longitudinal) [30] and the maximum change 
of 2.1% is taken as systematic uncertainty. The efficiency of the Q cut depends on the 
X(3872) decay model. The dipion mass spectrum obtained in this analysis does not 
have enough accuracy to discriminate between reasonable models. Comparing the results 
obtained with the X(3872) — )• J/ipp decay models used by CDF [9j and by Belle |7|, we 
evaluated a 1% systematic uncertainty on the Q-cut efficiency. 

Finally, differences in the trigger efficiency between data and simulation are studied 
using events triggered independently of the J/ip candidate. Based on these studies an 
uncertainty of 2.9% is assigned. 

6 Results and conclusion 

With an integrated luminosity of 34.7pb^^ collected by the LHCb experiment, the pro- 
duction of the X(3872) particle is observed in pp collisions at y/s = 7TeV. The product 
of the production cross-section and the branching ratio into J/ipTc'^TT' is 

a{pp ^ X(3872) + anything) i3(X(3872) ^ J/ipn+n-) = 5.4 ± 1.3 (stat) ± 0.8 (syst) nb , 

for X(3872) mesons produced (either promptly or from the decay of other particles) with 
a rapidity between 2.5 and 4.5 and a transverse momentum between 5 and 20 GeV/c. 

Predictions for the X(3872) — )■ J/tpii^Ti^ production at the LHC are available from 
a non-relativistic QCD model which assumes that the cross-section is dominated by the 
production of charm quark pairs with negligible relative momentum [13]. The calculations 
are normalized using extrapolations from measurements performed at the Tevatron. When 
restricted to the kinematic range of our measurement and summed over prompt production 
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and production from 6-hadron decays, the results of Ref. [T^ yield 13.0 ± 2.7nb, where 
the quoted uncertainty originates from the experimental inputs used in the calculation. 
This prediction exceeds our measurement by 2. 4a. 

After calibration using J/ip decays, the masses of both the X(3872) and 

il^i^S) mesons, reconstructed in the same J/iIjtt~^7t~ final state, are measured to be 

mx(3872) = 38 71 .95 ± 0.48 (stat) ±0.12 (syst) MeV/c^ , 
m^{2S) = 3686.12 ±0.06 (stat) ±0.10 (syst) MeV/c^ 

in agreement with the current world averages [23], and with the recent X(3872) mass 
measurement from Belle The measurements of the X(3872) mass are consistent, 
within uncertainties, with the sum of the and D*^ masses, 3871.79 ± 0.29 MeV/c^, 
computed from the results of the global PDG fit of the charm meson masses 
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